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Amongst iron regime nuclei, ,3-decay rates on titanium isotopes are considered to be 
important during the late phases of evolution of massive stars. The key /3-decay isotopes 
during presupernova evolution were searched from available literature and a microscopic 
calculation of the decay rates were performed using the proton- neutron quasiparticle ran- 
dom phase approximation (pn-QRPA) theory. As per earlier simulation results electron 
capture and /3-decay on certain isotopes of titanium are considered to be important for 
the presupernova evolution of massive stars. Earlier the stellar electron capture rates and 
neutrino energy loss rates due to relevant titanium isotopes were presented. In this paper 
we finally present the /3-decay rates of key titanium isotopes in stellar environment. The 
results are also compared against previous calculations. The pn-QRPA /3-decay rates 
are bigger at high stellar temperatures and smaller at high stellar densities compared 
to the large scale shell model results. This study can prove useful for the core-collapse 
simulators. 

PACS Number(s): 97.10.Cv, 26.50.+X, 23.40.-s, 23.40.Bw, 21.60.Jz 



1. Introduction 

Shortly after the discovery of the Gamow- Teller (GT) resonance, Bethe and collab- 
orators m suggested its importance for stellar weak- interaction mediated reactions. 
The calculations of stellar /3-decays and electron captures are very sensitive to the 
distribution and the total strength of the GT giant resonance. The centroid of the 
GT distribution functions determines the effective energy of the capture and decay 
reactions. Due to these weak interaction processes the value of Y e (lepton-to-baryon 
ratio) for a massive star changes from 1 (during hydrogen burning) to roughly 0.5 
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(at the beginning of carbon burning) and finally to around 0.42 just before the 
collapse to a supernova explosion. /3-decay are relatively unimportant at the be- 
ginning of carbon burning phase of massive stars (Y e ~ 0.5). However it becomes 
increasingly competitive for neutron-rich nuclei due to an increase in phase space 
related to larger values. In earlier simulations of core-collapse studies /3-decay 
of heavy nuclei were not considered. However it was soon acknowledged that during 
the early stages of the collapse, a strong /3-decay rate can contribute both to the 
cooling rate (via the stellar antineutrinos produced) and to a much larger value of 
Y e at the formation of homologous core. This in turn would result in a smaller enve- 
lope and to a more energetic shock wave. Aufderheide and collaborators ^ stressed 
on the importance of /3-decay rates in the iron core prior to the collapse. The au- 
thors also found in their calculation that the /3-decay rates were strong enough and 
at times larger than the competing electron capture rates for Y e ~ 0.42 — 0.46. 
The temporal variation of Y e within the core of a massive star has a pivotal role to 
play in the stellar evolution and a fine-tuning of this parameter at various stages 
of presupernova evolution is the key to generate an explosion. Collapse simulators 
world-wide find it challenging to transform the collapse of the core of massive stars 
to an explosion and to date there have been no successfully simulated spherically 
symmetric explosions. Even the 2D simulations (addition of convection) performed 
with a Boltzmann solver for the neutrino transport fails to convert the collapse into 
an ex plosio n (It is worth mentioning that recently a few simulation groups (e.g. 
Refs. I I I I) have reported successful explosions in 2D mode.) Additional energy 
sources (e.g. magnetic fields and rotations) were also sought that might transport 
energy to the mantle and lead to an explosion. World-wide core-collapse simulators 
are still working hard to come up with a convincing and decisive mode of producing 
explosions. 

The calculation of stellar weak interaction rates has been performed by a num- 
ber of authors during 1960's and 1970's. The rates were calculated using the best 
available physics of the time. However the first breakthrough was achieved in 1980 
when Fuller, Fowler and Newman (commonly referred to as FFN) ^ used a simple 
shell model to compute the location and strength of the GT resonances for 226 
nuclei in the mass range 21 < A < 60. FFN estimated the GT contributions to the 
rates by a parametrization based on the independent particle model. Aufderheide 
et al. ^ later extended the FFN work for heavier nuclei with A > 60 and took into 
consideration the quenching of the GT strength neglected by FFN. 

Experimentally it has been shown that for (p, n) and (n,p) reactions the 0° 
cross sections for such transitions are proportional to the squares of the matrix 
elements for the GT /3 decay between the same states (e.g. Ref. ^j). Results of 
such measurements disclosed that, in contrast to the independent particle model, 
the total GT strength is quenched and fragmented over many final states in the 
daughter nucleus caused by the residual nucleon-nucleon correlations. Both these 
effects are caused by the residual interaction among the valence nucleons and an 
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accurate description of these correlations is essential for a reliable evaluation of the 
stellar weak interaction rates due to the strong phase space energy dependence. 
The measured data from various (p,n) and (n,p) experiments also revealed the 
misplacement of the GT centroid adopted in the parameterizations of FFN and 
subsequently used in the calculation of weak rates by Ref. 1=1 Since then theoret- 
ical efforts were concentrated on the microscopic calculations of weak-interaction 
mediated rates of iron-regime nuclide. Two such widely used models are the large- 
scale shell model (LSSM)(e.g. Ref.®) and the prot on-n eutron quasiparticle random 
phase approximation (pn-QRPA) theory (e.g. Ref. EH). Both LSSM and pn-QRPA 
model later validated the finding of Ref. on microscopic grounds rather than phe- 
nomenological parametrization of GT centroids and strengths, that /3-decay rates 
are larger than electron capture rates for a certain range of Y e values. The LSSM 
calculation, however, makes use of the so-called Brink's hypothesis in the electron 
capture direction and back-resonances in the /3-decay direction to approximate the 
contributions from high-lying excited state GT strength distributions. Brink's hy- 
pothesis states that GT strength distribution on excited states is identical to that 
from ground state, shifted only by the excitation energy of the state. GT back res- 
onances are the states reached by the strong GT transitions in the inverse process 
(electron capture) built on ground and excited states. /3-decay rate calculation for 
iron-regime nuclei are also important during the late stages of stellar evolution. 

Isotopes of titanium are amongst the key iron-regime nuclei that play a key role 
in the developments during the late phases of stellar evolution. A search was per- 
formed by authors of Ref. 1=1 for most important electron capture and /3-decay nuclei 
after core silicon burning in massive stars. The authors assumed that the temper- 
atures within the iron core were high enough for the strong and electromagnetic 
reactions to be in equilibrium. The nuclear Saha equation was then used to com- 
pute the isotopic abundances. They tabulated the 71 top /3-decay nuclei averaged 
throughout the stellar trajectory for 0.40 < Y e < 0.5 (see Table 26 therein). From 
this table, six isotopes of titanium, namely 51 > 52 > 53 > 54 i 55 > 56 Ti, were short-listed as 
those whose /3-decay rates were suggested to have a significant contribution during 
and after core silicon burning phases in massive stars. A microscopic calculation of 
stellar /3-decay rates was performed for these titanium isotopes using the pn-QRPA 
theory. Because of the high temperatures prevailing during the presupernova and 
supernova phase of a massive star, there is a reasonable probability of occupation of 
parent excited states and the total weak interaction rates have a finite contribution 
form these excited states. The pn-QRPA theory allows a microscopic state-by-state 
calculation of all these partial rates and this feature of the model greatly enhances 
the reliability of the calculated rates in stellar matter. The pn-QRPA model can 
handle any arbitrarily heavy system of nucleons as it has access to a luxurious model 
space of up to Ifrw shells. The pn-QRPA model was successfull y used to calculat e 



weak interaction rates on important iron-regime nuclei (e.g. Refs J ' ' ' ' ' ). 
Earlier Nabi and collaborators^^ presented a detailed analysis of the calculation of 
stellar electron capture rates on twenty two titanium isotopes. Later a calculation 
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of neutrino and antineutrino energy loss rates in massive stars due to isotopes of 
titanium was also presented ^1 In this paper we present the /3-decay rates due to 
the above mentioned six isotopes of titanium in stellar matter. This completes the 
analysis of weak interaction rates of key titanium isotopes considered to be impor- 
tant for the presupernova evolution of massive stars. These microscopic rates can 
be of great utility for core-collapse simulators. The next section discusses briefly 
the formalism and presents the calculated /3-decay rates. Comparison with previous 
calculations is also presented in this section. Section 3 finally summarizes the main 
conclusions. 



2. Calculations and Results 



The Hamiltonian of the pn-QRPA model and its diagonalization were discussed 
earlier in Ref. The /3-decay rates of a transition from the i th state of the parent 
to the j th state of the daughter nucleus are given by 



\j - 



In 2 



D 



[B(F)ij + (gA/gvfBiGT)^} /g(T, p, E f ) 



(1) 



The value of D was taken to be 6295s . B(F) and B(GT) are reduced transition 
probabilities of the Fermi and Gamow- Teller (GT) transitions, respectively, 



B{GT\ 



2J 2 + 1 
1 

2Ji + l 



i > 



(2) 



< 



±z k IM>I 2 



(3) 



Here a k is the spin operator and t\ stands for the isospin raising and lowering 
operator. 

The ff- are the phase space integrals and are functions of stellar temperature 
(T), density (p) and Fermi energy (Ef) of the electrons. They are explicitly given 

by 



r 



wy w 2 



l{w m - w) d F(+Z,w)(l -G-)dw, 



(4) 



In Eq. (j4}, w is the total energy of the electron including its rest mass. w m is the 
total /3-decay energy, 



w.„ 



m d + Ei — Ea , 



(5) 



where m p and Ei are mass and excitation energies of the parent nucleus, and rrid 
and Ej of the daughter nucleus, respectively. F(+Z, w) are the Fermi functions and 
were calculated according to the procedure adopted by Gove and Martin 021. G_ 
is the Fermi-Dirac distribution function for electrons. Details of the calculation of 
reduced transition probabilities can be found in Ref. ^3. Construction of parent 
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and daughter excited states and calculation of transition amplitudes between these 
states can be seen in Ref. ^1 

The total /3-decay rate per unit time per nucleus is finally given by 



where Pi is the probability of occupation of parent excited states and follows the 
normal Boltzmann distribution. After the calculation of all partial rates for the 
transition i — > j the summation was carried out over all initial and final states until 
satisfactory convergence was achieved in the rate calculation. The pn-QRPA theory 
allows a microscopic state-by-state calculation of both sums present in Eq. ([6]). This 
feature of the pn-QRPA model greatly increases the reliability of the calculated 
rates over other models in stellar matter where there exists a finite probability of 
occupation of excited states. 

Stellar /3-decay rates are sensitive functions of the available phase space, (Qp + 
Ei — Ej). It is very much possible for nuclei with negative Qp values to undergo /3- 
decays in stellar environment (strictly forbidden under terrestrial conditions) as the 
phase space can become positive depending on the calculated energy eigenvalues of 
the underlying theoretical model. The calculated /3-decay rates of 51 > 52 . 53 . 54 . 55 > 56 Ti 
in stellar environment are presented in Table 1. The calculated rates are tabulated 
on an abbreviated density scale. The first column gives log(pY" e ) in units of gcm~ 3 , 
where p is the baryon density and Y e is the ratio of the electron number to the 
baryon number. Stellar temperatures (Tg) are stated in 10 9 K. The third column 
shows the calculated values of the Fermi energy in units of MeV. The last six 
columns give the calculated /3-decay rates of selected titanium isotopes in units 
of s~ x . The calculated /3-decay rates are tabulated in logarithmic (to base 10) 
scale. In the table, -100 means that the rate is smaller than 10 -100 s — ^ It can 
be seen from the table that at low stellar densities and temperatures, 56 Ti has 
the strongest /3-decay rate whereas 52 Ti the weakest. As the core stiffens from 
density pY e {gcm~ 3 ] = 10 to 10 4 the /3-decay rates does not change appreciably for 
a particular core temperature. As the core stiffens further the /3-decay rates start 
decreasing by orders of magnitude because of the considerable reduction in available 
phase space. The /3-decay rates increase monotonically with increasing temperature. 
Positron capture rates act in the same direction as /3-decay rates and at times tend 
to compete with the later. The positron capture rates on these titanium isotopes 
were also calculated. Table 2 shows the ratio of the calculated positron capture to 
/3-decay rates for isotopes of titanium at selected temperature and density points. It 
can be seen from this table that only at high stellar temperatures, Tg[-K] ~ 30, does 
the positron capture rate compete with the /3-decay rates. For all other physical 
conditions the positron capture rates are smaller by many orders of magnitude 
and hence can be safely neglected as compared to the /3-decay rates during the 
presupernova evolution of massive stars where the core temperature is considerably 
less. The complete electronic version (ASCII files) of these rates may be requested 
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Table 1. /3-decay rates of 51, 52, 53, 54,55,56^; f or se i ec t e d densities and temperatures in 
stellar matter. logpY" e has units of g/cm? , where p is the baryon density and Y e is the 
ratio of the lepton number to the baryon number. Temperatures (Tg) are given in units 
of 10 9 K. The calculated Fermi energy is denoted by Ef and is given in units of MeV. 
All calculated /3-dccay rates arc tabulated in logarithmic (to base 10) scale in units of 
s -1 . In the table, -100 means that the rate is smaller than lO -100 ^" 1 . 



logpY 


Tg 




51 T ; 


52 Ti 


53 Ti 


54 T; 


55 Ti 


56 Ti 


1.0 


0.01 


0.508 


-2.698 


-3.465 


-2.137 


-0.466 


-0.146 


1.004 


1.0 


0.10 


0.453 


-2.732 


-3.465 


-2.137 


-0.466 


-0.146 


1.004 


1.0 


0.20 


0.377 


-2.777 


-3.465 


-2. 135 


-0.466 


-0. 146 


1.004 


1.0 


0.40 


0.205 


-2.812 


-3.465 


-2.106 


-0.466 


-0.145 


1.004 


1.0 


0.70 


0.008 


-2.830 


-3.465 


-2.049 


-0.466 


-0.108 


1.004 


1.0 


1.00 


0.000 


-2.834 


-3.465 


-2.012 


-0.466 


-0.009 


1.004 


1.0 


1.50 


0.000 


-2.807 


-3.466 


-1.979 


-0.466 


0. 169 


1.004 


1.0 


2.00 


0.000 


-2.739 


-3.468 


-1.958 


-0.466 


0.292 


1.004 


1 


3 00 


000 


-2 557 


-3 482 


-1 798 


-0 468 


437 


1 003 


1.0 


5.00 


0.000 


-2.074 


-3.220 


-0.913 


-0.471 


0.594 


1.002 


1.0 


10.00 


0.000 


-1.031 


-0.946 


0.126 


0.139 


0.814 


1.269 


1.0 


30.00 


0.000 


0.221 


0.699 


1.172 


1.424 


1.616 


2.138 


4.0 


0.01 


0.523 


-2.699 


-3.468 


-2. 139 


-0.466 


-0. 146 


1.004 


4.0 


0.10 


0.516 


-2.733 


-3.467 


-2. 138 


-0.466 


-0. 146 


1.004 


4.0 


0.20 


0.498 


-2.778 


-3.466 


-2. 136 


-0.466 


-0. 146 


1.004 


4 


0.40 


0.444 


-2.813 


-3.466 


-2. 107 


-0.466 


-0. 145 


1.004 


4.0 


0.70 


0.337 


-2.830 


-3.466 


-2.050 


-0.466 


-0. 108 


1.004 


4.0 


1 00 


209 


-2 834 


-3 466 


-2 013 


-0 466 


-0 009 


1 004 


4.0 


1.50 


0.047 


-2.807 


-3.466 


-1.979 


-0.466 


0.169 


1.004 


4.0 


2.00 


0.014 


-2.739 


-3.468 


-1.958 


-0.466 


0.292 


1.004 


4.0 


3.00 


0.004 


-2.557 


-3.482 


-1.798 


-0.469 


0.437 


1.003 


4.0 


5.00 


0.001 


-2.074 


-3.220 


-0.913 


-0.471 


0.594 


1.002 


4.0 


10.00 


0.000 


-1.031 


-0.946 


0.126 


0.139 


0.814 


1.269 


4.0 


30.00 


0.000 


0.221 


0.699 


1.172 


1.424 


1.616 


2.138 


7.0 


0.01 


1.223 


-2.918 


-3.786 


-2.285 


-0.543 


-0. 192 


0.978 


7.0 


0.10 


1.222 


-2.946 


-3.786 


-2.285 


-0.543 


-0.192 


0.979 


7.0 


0.20 


1.222 


-2.986 


-3.785 


-2.283 


-0.543 


-0.192 


0.979 


7.0 


0.40 


1.219 


-3.015 


-3.783 


-2.242 


-0.543 


-0.191 


0.979 


7.0 


0.70 


1.212 


-3.027 


-3.777 


-2.165 


-0.542 


-0.151 


0.979 


7.0 


1.00 


1.200 


-3.024 


-3.768 


-2.116 


-0.541 


-0.044 


0.979 


7.0 


1.50 


1.173 


-2.980 


-3.749 


-2.069 


-0.538 


0.143 


0.980 


7.0 


2.00 


1.133 


-2.886 


-3.726 


-2.037 


-0.535 


0.272 


0.981 


7.0 


3.00 


1.021 


-2.658 


-3.681 


-1.845 


-0.528 


0.422 


0.983 


7.0 


5.00 


0.698 


-2.109 


-3.275 


-0.924 


-0.507 


0.585 


0.989 


7.0 


10.00 


0.196 


-1.037 


-0.950 


0.122 


0.134 


0.811 


1.266 


7.0 


30.00 


0.021 


0.221 


0.699 


1.171 


1.423 


1.615 


2.138 


10.0 


0.01 


11.118 


-100 


-100 


-100 


-100 


-100 


-100 


10.0 


0.10 


11.118 


-100 


-100 


-100 


-100 


-100 


-100 


10.0 


0.20 


11.118 


-100 


-100 


-100 


-100 


-87.308 


-93.598 


10.0 


0.40 


11.118 


-100 


-100 


-74.497 


-85.977 


-45.263 


-48.401 


10.0 


0.70 


11.117 


-62.404 


-65.400 


-43.640 


-50.168 


-26.906 


-28.699 


10.0 


1.00 


11.116 


-44.379 


-46.217 


-31.133 


-35.642 


-19.400 


-20.626 


10.0 


1.50 


11.113 


-30.186 


-31.101 


-21.236 


-24.133 


-13.413 


-14.125 


10.0 


2.00 


11.110 


-22.978 


-23.418 


-16.171 


-18.229 


-10.322 


-10.715 


10.0 


3.00 


11.099 


-15.627 


-15.575 


-10.948 


-12.114 


-7.096 


-7.084 


10.0 


5.00 


11.063 


-9.553 


-9.075 


-6.539 


-6.902 


-4.308 


-3.841 


10.0 


10.00 


10.898 


-4.718 


-3.878 


-2.894 


-2.530 


-1.871 


-0.907 


10.0 


30.00 


9.163 


-0.701 


-0.090 


0.362 


0.693 


0.866 


1.471 


11.0 


0.01 


23.934 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


0.10 


23.934 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


0.20 


23.934 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


0.40 


23.934 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


0.70 


23.934 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


1.00 


23.933 


-100 


-100 


-95.730 


-100 


-83.987 


-85.027 


11.0 


1.50 


23.932 


-73.255 


-74.170 


-64.305 


-67.176 


-56.474 


-57.036 


11.0 


2.00 


23.930 


-55.285 


-55.725 


-48.478 


-50.515 


-42.622 


-42.893 


11.0 


3.00 


23.925 


-37.175 


-37.122 


-32.496 


-33.647 


-28.639 


-28.536 


11.0 


5.00 


23.908 


-22.501 


-22.023 


-19.487 


-19.840 


-17.252 


-16.709 


11.0 


10.00 


23.832 


-11.236 


-10.393 


-9.411 


-9.033 


-8.380 


-7.344 


11.0 


30.00 


23.016 


-2.980 


-2.334 


-1.889 


-1.527 


-1.361 


-0.703 
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Table 2. Ratio of calculated positron capture rates to /3-decay rates for 51, 52, 53, 54 > 55 " 56 Ti for selected 
densities and temperatures in stellar matter. logpY e has units of g/cm' i , where p is the baryon density 
and Y e is the ratio of the lepton number to the baryon number. Temperatures (Tg) are given in units 
of 10 9 K. The calculated Fermi energy is denoted by Ef and is given in units of MeV. Cases with no 
reported ratio imply that either (or both) of the calculated rate(s) is (are) smaller than 10 — 100 s _1 . 



logpY,, 


Tg 


Ef 


51 Ti 


52 T; 


53 Ti 


54 T ; 


55 Tj 


66 Ti 


1.0 


0.10 


0.453 


1.20E-53 


5.42E-52 


3.18E-53 


5.82E-54 


2.89E-54 


1.04E-54 


1.0 


1.00 


0.000 


2.70E-05 


1.78E-03 


7.21E-05 


1.50E-05 


5.52E-06 


2.51E-06 


1.0 


3.00 


0.000 


1.46E-02 


2.32E+00 


5.20E-02 


1.23E-02 


1.79E-03 


1.75E-03 


1.0 


10.00 


0.000 


1.25E+00 


2.49E+00 


4.20E-01 


7.85E-01 


1.94E-01 


1.69E-01 


1.0 


30.00 


0.000 


2.25E+02 


5.55E+01 


6.01E+01 


3.07E+01 


3.02E+01 


1.15E+01 


4.0 


0.10 


0.516 


8.04E-57 


3.64E-55 


2.13E-56 


3.90E-57 


1.93E-57 


6.98E-58 


4.0 


1.00 


0.209 


2.41E-06 


1.59E-04 


6.44E-06 


1.34E-06 


4.92E-07 


2.24E-07 


4.0 


3.00 


0.004 


1.44E-02 


2.29E+00 


5.13E-02 


1.21E-02 


1.76E-03 


1.73E-03 


4.0 


10.00 


0.000 


1.25E+00 


2.49E+00 


4.20E-01 


7.85E-01 


1.94E-01 


1.69E-01 


4.0 


30.00 


0.000 


2.25E+02 


5.55E+01 


6.01E+01 


3.08E+01 


3.02E+01 


1.15E+01 


7.0 


0.10 


1.222 


3.37E-92 


1.95E-90 


7.67E-92 


1.19E-92 


5.51E-93 


1.90E-93 


7.0 


1.00 


1.200 


3.75E-11 


3.21E-09 


8.20E-11 


1.60E-11 


5.37E-12 


2.38E-12 


7.0 


3.00 


1.021 


3.62E-04 


7.18E-02 


1.13E-03 


2.77E-04 


3.64E-05 


3.61E-05 


7.0 


10.00 
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from the corresponding author. 

The calculation of /3-decay rates was also compared with previous calculations. 
For the sake of comparison we considered the pioneer calculations of FFN El and 
those performed using the large-scale shell model (LSSM) Figure [1] depicts the 
comparison of /3-decay rates of 51 Ti with earlier calculations. The upper panel dis- 
plays the ratio of calculated rates to the LSSM rates, R/3(QRPA/LSSM), while 
the lower panel shows the corresponding comparison between FFN and LSSM cal- 
culations, R[j{FFN/ LSSM). All graphs are drawn at four selected values of stellar 
densities (pY e [gcm~ 3 ] = 10 1 ,10 4 ,10 7 and 10 11 ). These values correspond roughly 
to low, medium-low, medium-high and high stellar densities, respectively. The se- 
lected values for temperature on the abscissa are Tg[A'] = 1, 1.5, 2, 3, 5, 10 and 30. 
It can be seen from Figure[T]that the pn-QRPA calculated /3-decay rates of 51 Ti are 
generally in good agreement with LSSM rates (within a factor 5). The comparison 
between FFN and LSSM calculations is also good except at high temperatures and 
densities. FFN did not take into effect the process of particle emission from excited 
states and their parent excitation energies extended well beyond the particle decay 
channel. The occupation probability of these high lying excited states becomes fi- 
nite at high temperatures and densities. FFN rates are consequently bigger by more 
than one order of magnitude at Tg[A] = 30. 

For the case of /3-decay rates of 52 Ti, the LSSM rates are slightly bigger (Fig- 
ure[2]) compared to the reported decay rates. At low densities the LSSM decay rates 
are bigger by up to a factor of five whereas at high densities they are bigger by as 
much as factor of twenty. Whereas the individual discrete transitions between ini- 
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Fig. 1. (Color online) Ratio of reported /3-decay rates to those calculated using LSSM (upper 
panel) for 5 Ti as function of stellar temperatures and densities. Tg gives the stellar temperature 
in units of 10 9 K. In the legend, logpY e gives the log to base 10 of stellar density in units of 
gcm~ 3 . The lower panel shows the corresponding ratio of FFN /3-decay rates to those calculated 
using LSSM. 



tial and final states matter at low temperatures and densities, it is the total B(GT) 
strength that matters at high temperatures and densities. Brink's hypothesis is not 
assumed in this calculation (which was adopted in LSSM calculation of weak rates). 
At high temperatures, Tg[l£] = 30, the agreement is very good hinting towards the 
fact that the total B(GT) strength calculated in both models match reasonably 
well. The lower panel shows that FFN rates are again in reasonable comparison 
with LSSM rates except at high temperatures for reasons mentioned above. 

The agreement with LSSM /3-decay rates for the odd-A nucleus, 53 Ti, is reason- 
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Fig. 2. (Color online) Same as Fig.l. but for /3-decay of 52 Ti 

able (within a factor 5) and depicted in Figure |3] However the pn-QRPA rates are 
bigger by a factor of thirteen at high temperatures. Comparison between FFN and 
LSSM calculations (lower panel) reveals that FFN /3-decay rates are much bigger 
specially at high temperatures. Authors in Ref. ^ reported that, for even-even and 
odd-A nuclei, FFN systematically placed the back resonance at much lower ener- 
gies and concluded that contribution of the back resonances to the /3-decay rates 
for these nuclei decreases. They estimated that LSSM /3-decay rates as a result were 
smaller, on the average, by a factor of 20 (40) as compared to the FFN /3-decay 
rates for even-even (odd-A) nuclei. 

FFN did not calculate the /3-decay rates of 54 Ti and as such a mutual compari- 
son of the three calculations was not possible for this even-even isotope of titanium. 
We note from Figure |4] that our /3-decay rates of 55 Ti are slightly bigger (factor of 
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Fig. 3. (Color online) Same as Fig.l. but for /3-decay of 53 Ti 



2 - 3) than the LSSM rates at pY e [gcm- 3 } = 10\ 10 4 ,10 7 . At high densities the 
LSSM rates exceed the pn-QRPA rates as in previous cases. At high stellar temper- 
atures, the reported rates are more than an order of magnitude bigger than LSSM 
/3-decay rates. At higher temperatures excited state GT strength distributions are 
required for the calculation of weak rates (parent excited states have a finite prob- 
ability of occupation). The LSSM employed the Brink's hypothesis in the electron 
capture direction and back-resonances in the /3-decay direction to approximate the 
contributions from high-lying excited state GT strength distributions. On the other 
hand the pn-QRPA model performs a microscopic calculation of the GT strength 
distributions for all parent excited states and provides a fairly reliable estimate of 
the total stellar rates. The lower panel of Figure |4] shows a whopping seven orders 
of magnitude bigger FFN decay rates at low temperatures and high densities. How- 
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Fig. 4. (Color online) Same as Fig.l. but for /3-decay of 55 Ti 



ever for the same temperature and density domain the two microscopic calculations 
(pn-QRPA and LSSM) are in good agreement hinting towards the fact that FFN 
overestimated their /3-decay rates by orders of magnitude. The centroids of the GT 
distribution functions determine the effective energy of the capture and decay reac- 
tions. FFN estimated the GT contributions to the rates by a parametrization based 
on the independent particle model and estimated the GT centroids using zeroth- 
order (0hu> ) shell model. Stellar decay rates are fragile functions of the available 
phase space, (Qp + Ei — Ej). It is worth mentioning that these /3-decay rates can 
change by orders of magnitude by a mere change of 0.5 MeV, or less in available 
phase space and are more reflective of the uncertainties in the calculation of en- 
ergies. It is to be noted that whereas pn-QRPA and LSSM used the value of 7.34 
MeV and 7.81 MeV, respectively, as Q-value for this decay reaction, FFN used a 
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much larger value of 8.64 MeV for the same reaction. 

The comparison of LSSM and pn-QRPA /3-decay rates of 56 Ti is similar to the 
previous case of 55 Ti (see Figure [5]). At Tgfif] = 30 the reported decay rates are 




Fig. 5. (Color online) Same as Fig.l. but for /3-decay of 56 Ti 



around a factor of thirty bigger than LSSM rates. For a change we notice in the 
lower panel that the FFN decay rates of 56 Ti are much smaller than the LSSM 
rates. The rates are more than six orders of magnitude smaller at low temperatures 
and densities. For the same physical conditions the reported rates are again in far 
better agreement with the LSSM numbers again hinting toward some problems in 
the FFN calculations. It was found by authors in Ref. ^that the approach of FFN 
was not always reliable in its estimates of the location of GT strength. Unmeasured 
matrix elements for allowed transitions were assigned an average value of log ft =5 
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in FFN calculations. Measurements were performed for the /3-decay of 56 Ti 
after the FFN calculations. This measurement was properly incorporated in the 
microscopic calculation of pn-QRPA (this work) and LSSM and resulted in a much 
bigger /3-decay of 56 Ti. 



3. Conclusions 

Titanium isotopes are amongst the key iron-regime nuclei that play a significant 
role in the late phases of stellar evolution of massive stars. The /3-decay rates of 
these isotopes increase the value of the lepton-to-baryon fraction (Y e ) during the 
late phases of stellar evolution. The temporal variation of Y e within the core of a 
massive star has a pivotal role to play in the stellar evolution and a fine-tuning of 
this parameter at various stages of presupernova evolution is the key to generate an 
explosion. Six isotopes of titanium, 51 ^ 52 > 53 < 54 > 55 < 56 Ti, were short-listed as important 
/3-decay nuclei as the product of their abundance and /3-decay rates can cause a 
substantial change in the time rate of change of lepton-to-baryon fraction. The pn- 
QRPA model was used to calculate the /3-decay rates of these six titanium isotopes 
in a microscopic fashion. The pn-QRPA model has two distinct and important 
advantages as compared to other models. Firstly it can handle any arbitrarily heavy 
system of nucleons since the calculation is performed in a luxurious model space 
of up to 7 major oscillator shells. Further it is the only available model that can 
calculate all excited state GT strength distributions in a microscopic fashion which 
greatly increases its utility in stellar calculations. 

The /3-decay rates were calculated on a detailed density-temperature grid point. 
The positron capture rates on titanium isotopes were also calculated and it was 
shown that during the presupernova evolution of massive stars the positron capture 
rates can be safely neglected in comparison with the /3-decay rates. The ASCII 
files of the calculated rates can be requested from the corresponding author. The 
pn-QRPA /3-decay rates were also compared against previous calculations. Our 
study validates the finding of authors of Ref. that FFN systematically placed 
the back resonance at much lower energies. Consequently the FFN overestimated 
the /3-decay rates for even-even and odd-A nuclei. The comparison with LSSM 
calculation is generally fair. However at high densities the LSSM rates are enhanced 
whereas at high temperatures the pn-QRPA /3-decay rates are much bigger. The 
study also stresses on the fact that the Brink's hypothesis and back resonances 
(employed in previous calculations) are not a good approximation and it is desirable 
to microscopically calculate all excited state GT strength distribution functions for 
a reliable estimate of stellar /3-decay rates. Core-collapse simulators might find it 
useful to employ the pn-QRPA /3-decay rates as an alternate (microscopic and 
reliable) nuclear physics input in their codes. 
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